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ABSTRACT 

We report molecular line observations in the CO J=3— >2, 6—7-5 and 7— >6 
transitions, made using the Atacama Pathfinder Experiment Telescope (APEX), 
toward the massive and dense core IRAS 16562—3959. This core harbors a string 
of radio sources thought to be powered by a central collimated jet of ionized gas. 
The molecular observations show the presence of high velocity gas exhibiting a 
quadrupolar morphology, most likely produced by the presence of two collimated 
outflows. The southeast-northwest molecular outflow is aligned with the string 
of radio continuum sources, suggesting it is driven by the jet. We find that the 
excitation temperature of the gas in the SE-NW outflow is high, with values of 145 
and 120 K for the blueshifted and redshifted lobes, respectively. This outflow has 
a total mass of 1.92 M , a total momentum of ~ 89 M km s -1 and an averaged 
momentum rate of ~ 3.0 x 10~ 2 M Q km s _1 yr _1 , values characteristics of flows 
driven by young massive stellar objects with high luminosities (L ho \ ~ 2 x 10 4 
L ). Complementary data taken with the Atacama Submillimeter Telescope 
Experiment (ASTE) in high density and shock tracers support the picture that 
IRAS 16562—3959 is an accreting young massive star associated with an ionized 
jet, which is the energy source of a molecular outflow. 

Subject headings: ISM: individual (IRAS 16562-3959) — ISM: jets and outflows - 
radio continuum: stars — stars: formation 
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INTRODUCTION 



IRAS 16562-3959 (also G345.4938+01.4677) is a mid-infrared (MIR) source, with 
MSX (Mid course Space Experim ent) colors characteristics of massive young stellar objects 



(MYSOs) (ILumsden et al. 



20021 ). Observations of dust continuum emission at 1.2 mm, 



made with SIMBA at the SEST telescope, show that IRAS 16562—3 959 is associated w ith 



a massive (1.3 x 1O 3 M ) and dense (9 x 10 5 cm 3 ) mo lecular core (Lopez et al.ll201ll ). 



Maser emission has been detected in transitions of OH rtCaswelllll998 



20041 1. giving further 



support to the hypothesis that IRAS 16562—3959 is a young high-mass star forming region. 



Recently, iGuzman et al.l (120101 ) reported the detection toward IRAS 16562—3959 of a 
string of radio continuum emission consisting of a bright compact central object and four 
outer lobes. They argue that the continuum emission from the central object corresponds to 
free-free emission from a thermal jet whereas the radio emission from the lobes correspond 
to thermal emission arising in shocks generated from the interaction of the collimated wind 
with the surrounding medium. Ass uming that IRAS 165 62—3959 is located at a distance of 



1.6 ± 0.1 kpc (Vlsr = -12.6 km/s, 



Urquhart et al 



2008 



Casoli et al 



;he to tal far-infrared luminosity, 
19861 ). is ~ 7.0 x 10 4 L . This 



computed using the observed IRAS fluxes (see 
implies that IRAS 16562—3959 harbors the most luminous MYSO known to date associated 
with an ionized jet. 

Mid-infrared imag es taken from the Spitzer Space Telescope-GLIMPSE survey data 



( [Benjamin et al 



2003) show bright emission associated with the jet in all IRAC bands. In 
particular, at 8.0 /im there is extended emission roughly aligned with the jet axis. Further, 
there appears t o be an excess of emission at 4.5/xm , which is thought to be indicative of 



shock activity ( jCyganowski et al. 



2008 



Chambers et al. 



2009|). The GLIMPSE data also 



x The near/far kinematic distance ambiguity was resolved in 



Faundez et al. 



(120041 ) 
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show that there is at least one more young stellar object (YSO) within the core, located 
about 15" northeast of the jet, that is not detected at 2 /im and hence likely to be deeply 
embedded. 

In this paper we present molecular line observations toward IRAS 16562—3959 in 
the CO (3 — > 2), CO (6 — > 5), and CO (7 — > 6) transitions, made using APEX, which have 
revealed the presence of two high velocity bipolar molecular outflows. Their observed and 
derived characteristics suggest that the driving sources are young, luminous protostellar 
objects. One of the outflows is aligned in the same direction as the string of radio sources 
strongly suggesting that it is being driven by the thermal radio jet. 



2. OBSERVATIONS 



The observations were made using the 12-m Atacama Pathfinder Experiment Telescope 
(APEX) located at Llano de Chajnantor, Chile. A detailed description of APEX and 



its performance are given by 



Glisten et al. 



(120061 ). The observed transitions and basic 



observational parameters are summarized in Table [TJ Columns 1 and 2 give, respectively, 
the observed transition and line frequency. Columns 3, 4 and 5 give the telescope used, the 
half-power beam width and main beam efficiency at the observed frequency. Columns 6 to 9 
give, respectively, the number of positions observed, the angular spacing, the channel width, 
and the resulting rms noise in antenna temperature, for each of the observed transitions. 

The CO (3 — > 2) observations were made during June, 2008 and July, 2009. The 



front-end consisted o: 



band flRisacher et al. 



Spectrometer (IKlein et al.l 



a sin gle pixel heterodyne SiS receiver operating in the 275-370 GHz 



2006 



For the back-end we used the APEX Fast Fourier Transform 



20061 ) with a bandwidth of 1 GHz and 2048 frequency channels. 



The velocity coverage was ~ 870 km s 1 centered at —13.6 km s 1 Vlsr with a channel 
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width of 0.42 km s -1 . In this transition the emission was mapped within a region of 
100" x 100" in size, with 20" angular spacing (roughly full beam spacing), centered on 
the position of the radio source. System temperatures were typically 270 K. The pointing 
accuracy is 2 - 3" calibrated towards the source RT-Sco {a = 17 h 03 m 32.6 s , 5 = -36 °55'14") 
and the uncertainty in the absolute flux scale is ~ 10%. The observations were performed in 
position-switching mode under good atmospheric conditions (r ~ 0.07) using a = 16 /l 58 m 15 s 
and 5 = —39 °00'00"as off-position. The integration time on-source in each position was 
~ 30 seconds, resulting in an rms noise of typically 0.1 K in antenna temperature per 
channel. 

The CO(6 —> 5) and CO(7 — > 6) observations were made during August, 2009, under 
good weather conditions (0.35-0.45 mm of water vapor, corresponding to 0.5 < r < 0.8) 
using an off-position located at Aa = +6 m . The pointing calibration source was NGC6334, 
attaining an accuracy of ~ 2". The front-end consisted of a dual heterodyne SIS receiver 
array of 7 pixels operating in the 600-720 GHz and 750-950 GHz atmospheric windows, 



known as the CHAMP + receiver rtGiisten et al. 



20081 ). In these two transitions we mapped 



the emission within a region of 110" x 90" in size, with 4" angular spacing. The back-ends 
were operated with 2.4 GHz total bandwidth with 2048 channels each. 

We complement the above data with additional spectra taken toward the radio source 
position in the SiO(8^7), HCO+(4^3) and C 18 0(3^2) lines. The SiO(8^7) observations 
were performed using APEX with the same APEX FFT spectrometer. The integration 
time was ~ 450 seconds giving an rms noise of 0.04 K. 



The HCO + (4— >3) and C 18 0(3— >2) spectra were taken as part of a 



arger 



mapping 



200J) during 



of the core using the ASTE 10 m submillimeter telescope (lEzawa et al. 
June- July, 2010, under good atmospheric conditions (r ~ 0.08). The front-end consisted in 
a heterodyne single pixel SiS receiver operating between 324 and 372 GHz. As back-end we 
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used the MAC XF-type digital spectro-correlator with 1024 channels. For the C 18 0(3— >2) 
line we used a total bandwidth of 512 MHz, giving a velocity resolution of about 0.43 
km s _1 , whereas for the HCO + (4— >3) transition we used a 125 MHz bandwidth, giving a 
velocity resolution of ~ 0.11km s _1 . A detailed description of these observations will be 
given elsewhere. 

3. RESULTS 

Figure [1] shows a grid of the CO spectra in the three observed transitions within a 
region of 100" x 100" in size, with 20" spacing. The spectra of the CO(6— )-5) and CO(7— >Q) 
transitions were convolved to match the angular resolution of the CO (3— >2) spectra. Clearly 
seen in this figure is the presence of broad and strong wing emission across the molecular 
core. 

Figure [2] shows the HCO+(4^3), C 18 0(3^2) and SiO(8^7) spectra observed toward 
the peak position of IRAS 16562—3959. The HCO + (4— >3) spectrum shows a double-peaked 
line profile, with a bright blue-shifted peak at —14.0 ± 0.2 km s _1 and a weaker red-shifted 
peak at —10.1 ± 0.2 km s _1 . On the other hand, the C 18 0(3— >2) line profile shows a 
symmetric single component with a peak line center velocity of —12.5 km s _1 . We also note 
that the HCO + (4— >3) and the SiO(8— >7) spectra show prominent high velocity wings. 

The 12 CO profiles exhibit redshifted wing emission mainly towards the west and south, 
up to an LSR velocity of 14 km s _1 , whereas blueshifted wing emission is seen mainly 
towards the east and north, up to an LSR velocity of —46 km s _1 . The full velocity range 
of the wing emission observed is then ~ 60 km s . The radial flow velocity is defined as 
I visr — v am b | , where v am b is the systemic velocity of the ambient gas. For the later we 
adopt a value of —12.5 km s -1 , corresponding to the peak velocity of the C 18 0(3— ^2) line 
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profile observed toward the central position of the core. The maximum outflow velocity is 
then ~ 33.5 km s _1 towards the blue and ~ 26.5 km s" 1 towards the red. Figure [3] shows 
contour maps of the wing emission in the three CO lines integrated over the velocity range 
—46 < vi sr < —22.0 km s _1 (blue contours) and — 3 < v;, st . < 14 km s -1 (red contours). 



The star marks the position of the radio jet ( iGuzman et al. 



20101 ). These velocity limits 



define the ranges in which emission from the high-velocity outflow gas is seen. We note that 
the blueshifted flow emission is slightly stronger and reaches higher radial flow velocities 
compared to the redshifted emission, hence the wider blueshifted velocity range. The inner 
limits of the ranges were selected as to leave out the contribution of the ambient cloud and 
chosen symmetrically respect to the ambient cloud velocity. 

The spatial distribution of the blueshifted and redshifted emissions shows a quadrupolar 
morphology, with blueshifted emission seen mainly toward the east and north and redshifted 
emission seen toward the west and south. We propose that the quadrupolar morphology 
is produced by the superposition of two bipolar outflows, one along a southeast-northwest 
direction (hereafter the SE-NW flow) with a position angle of ~ 107°, and a second outflow 
aligned roughly in a north-south (N-S) direction. The collimation factors — length divided 
by width — of the molecular outflows are not high: for the blueshifted and redshifted lobes 
of the SE-NW outflow are > 2 and ~ 1, respectively, whereas for the lobes of the N-S 
outflow we found values close to 1. These values are similar to those derived for massive 



outflows (IBeuther et al. 



20021 ) . We note, however, that the highest velocity gas emission of 



the SE-NW outflow appears to be associated with more collimated structures. 

Figure H] presents a contour map of the outflow emission in the CO (6— >5) line 



overlayed with 


the 8. 


(Guzman et al. 


2010) 



20101 ). It appears that the SE-NW outflow is associated with the string of 



radio sources, the peak position of the blue-shifted and red-shifted lobes being symmetrically 
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displaced from the bright central radio source. The SE blueshifted and NW redshifted lobes 
extend up to ~ 27" and ~ 32" from the central radio source, respectively. The symmetry 
axis of the SE-N W outflow is along a direction with a position angle of ~ 107°, roughly 



the sa me as the P.A. of the symmetry axis of the jet of 110°. As noted in 



Guzman et al. 



(120 10l ). the radio lobes and the jet are not completely aligned, showing a small bending, 



which is also se en in the SE-NW outflow. Possible bending mechanisms of protostellar jets 



are discussed in 



Fendt fc Zinneckerl (119981 ). but with the available data we can not discern 



between the various alternatives. 

The spectroscopic signatures of the HCO + (4— )-3) and C 18 0(3— >2) transitions suggest 
that the bulk of the molecular gas toward IRAS 16562—3959 is undergoing large-scale 



inward motions (e.g. 



Sanhueza et al. 



3. 



20101 ) . Infalling motions traced by optically thick 



molecular lines are expected to produce line profiles showing bl ue asymmetry, where as 



optically thin lines are expected to exhibit symmetrical profiles (IMardones et al 



1997J) 



ANALYSIS AND DISCUSSION 



4.1. K s -band evidence of shocked H 2 



Figure [5] presents an image of the Two Micron All Sky Survey (2MASS. ISkrutskie et al. 



20061 ) A^-band emission across an 8' x 8' region of the sky, centered near IRAS 16562—3959. 
Clearly seen toward the center is diffuse emission along the SE-NW direction extending by 
more than 1.5' on each side of the bright radio source. The position angle of this diffuse 
iC,-band emission is ~ 104°, a value similar to that of the SE-NW molecular outflow. We 
suggest that the outer parts of this diffuse emission, which are at greater distances from the 
central jet than the molecular lobes, are tracing shocked H 2 -2.12yum gas produced by an 
older episode of mass ejection from the jet. 



We note that the NIR (near infrared) emission detected close to the central jet comes 
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mainly from a region associated with the SE blueshifted lobe. Much less emission is seen 
associated with the NW redshifted lobe (see also Fig. 4 from Guzman et al. 2010). A 
similar situati on has been observed in other bipolar outflows, for example in BHR 71 



( jBourke 



20021 ) . Most likely, the strong NIR emission from the eastern, closer to the jet, 



feature correspond to scattered light from the inner walls of the outflow cavity. 

The second outflow detected toward IRAS 16562—3959 is roughly aligned in the 
north-south direction, having a symmetry axis with a position angle of 7°. The blueshifted 
lobe extends up to ~ 21" to the north of the jet and the redshifted lobe up to ~ 17" to the 
south. Interestingly, Fig. [5] also shows the presence of diffuse i^ s -band emission roughly 
along the N-S direction but located about 3.4' north and 2' south from the NS molecular 
outflow. This emission exhibits a bow shock morphology, characteristic of HH objects, 
most likely produced by the interac t ion o f a collimated flow with the ambient medium. 



However, as noted in 



Guzman et al 



(120 lOl ). this diffuse emission is bright in all mid- infrared 
MSX bands. In consequence, it is likely that this i^-band emission has an important NIR 
continuum component, in addition to the H2-2.12/zm line. 

We suggest that there is an independent source of energy driving th e NS outflow, 



althou gh from the available data we are not able to pinpoint its location. iGuzman et al. 



(]2010l ) have already argued, on the basis that a re-orientation of the central jet in ~ 90 



seems physically unfeasible, for the presence of two high-mass YSOs driving the outflows. 



4.2. Parameters of the outflows 



To compute physi cal parameters of th e molecular out f 



forma lism described in iBourke et al. 



(1997) 



Garden et al. 



ows we fo 



fll99lh and 



lowed the standard 



Goldsmith fc Langer 



(119991 ). assuming that the high velocity gas is optically thin and its excitational state can be 
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described by a single excitation temperature. If the transitions are sub-thermally excited, 
then the derived excitation temperature would correspond to a lower limit of the kinetic 
temperature of the outflowing gas. 



A gen eral discussion of the sou rces of errors have been given by 



Margulis fc Lada 



(119851 ) and ICabrit fc Bertoutl (119901 ) . The main sources of error arise from the difficulty 
in determining the contribution to the outflow in the velocity range of the ambient cloud, 
and not knowing the flow inclination. To be conservative we adopt as velocity boundary 
between the blue and red wing emission and the ambient emission the values of —22.0 and 
—3.0 km s _1 , respectively. 



4-2.1. Column densities and excitation temperatures 



The column density of CO molecules in the velocity range [t>l,t>2] is given by 

M „ 1v in i4 (T ex + hB/3k) ex P (Ej/kT cx ) T 2 
iV C o = 2.31 x 10 — - — ■ — r / r v dv cm , 



1 - exp (-hu/kT ex ) (J + 1) J vl 

where T ex is the rotational excitation temperature in K, J is the rotational quantum 
number of the lower state, Ej = hBJ(J + 1) is the energy of level J, where B is the 
rotational constant of the CO molecule. The frequency of the CO( J+l— » J) transition 
is v fa 2B(J + 1), and t v is the opacity of the material moving at velocity v (the latter 
measured in km s _1 ). Assuming a beam filling factor of 1, the observed main-beam 
brightness temperature is related to the opacity by 

r*t» = [M",T«) - J^,T bg )](l - e--) , (2) 

where J r (u,T) is defined as 

J r (v,T) - 



exp(hv/kT) — 1 
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In the optically thin limit, 

T* b (v) = [J r (v,T ex )-J r (v,T hg )]T v . (3) 

From the observations of the emission in the three transitions of CO it is possible to 
determine both the excitation temperature and column density of the outflowing gas using 
Eqs. ([1]) and and taking Tb g = 2.7 K. Figure [6] presents plots of the velocity integrated 
emission versus the J quantum number from the blueshifted and redshifted lobes of both 
outflows and from the ambient gas. The velocity intervals of integration for the blueshifted 
and redshifted emission are [—46, —22] and [—3, 14] km s _1 , respectively. The areas over 
which the emission has been spatially integrated for the different outflow lobes are shown 
in Fig. HJ In the positions in which lobes overlap we assumed equal contributions from 
each lobe, except at offset position (Act = +20" AS = +20") where we assumed that the 
contribution to the blueshifted emission from the NS and SE-NW outflows are in a 2:1 
ratio. Dashed lines show the results of the best fit to the observed integrated emission 
assuming optically thin conditions and filling factors of 1. The derived CO column densities 
and excitation temperatures are given in Table [2j The parameters for the ambient cloud 
were derived using the spectra observed at offset position (Aa = +40", AS = —40"), which 
does not show evidence of self-absorption nor high velocity gas. 

One of the less well known physical parameter of the gas in molecular outflows is its 
kinetic temperature. Our observations of three lines of CO allowed us to determine the 
excitation temperature of both outflows, concluding that their kinetic temperatures are 
high. The derived excitation temperatures of the high- velocity gas are 145 and 120 K for the 
blueshifted and redshifted lobes of the SE-NW outflow, and 76 and 88 K for the blueshifted 
and redshifted lobes of the NS outflow. All these temperatures are considerably higher than 
the ambient cloud temperature, indicating that part of the energy used in accelerating the 
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?as has also heated it. 



4-2.2. Mass, momentum and momentum rates 



The mass in the outflows can be computed from the derived column densities as, 



M=[H 2 / 12 CO] fi m J2 N co dA 



(4) 



where N C o is given by Eq. ([T]), fi m is the mean molecular mass per H 2 molecule, [H 2 / 12 CO] 
the molecular hydrogen to carbon monoxide abundance ratio, dA the size of the emitting 
area in an individual position, and the sum is over all the observed positions. 

The derived masses of the SE-NW and NS flows are given in Columns 2-4 of Table |3j 
Column 2 gives the mass of the outflows in the high-velocity (HV) ranges, computed using 
the column densities given in Tabled fi m = 2.3 m H , and [H 2 / 12 CO] = 10 4 . To estimate 
the contribution to the mass from the outflowing gas emitting in the same velocity range 
as the ambient cloud gas, referred as t he low- velocity (LV) outflow gas, we followed th e 



prescription of 



Margulis fc Ladal (119851 ). Using expression (A16) of 



Bourke et al. 



(119971 ) 



adopting as velocity boundary between the wing and ambient emissions values of —22 
km s _1 in the blue side and —3 km s _1 in the red side, we estimate that the mass of the 
low velocity flow is 1.26 and 0.92 M for the SE-NW and NS flow, respectively. The total 
masses for the SE-NW and NS outflows are then 1.92 and 1.32 M , respectively. 

The momentum of the gas can be estimated from the first moment of the line emission. 
This procedure has the advantage of using the detailed information of the spectra. The 
derivation of flow kinematic an d dynamical param eters by using moments of the spectrum 



has been described in detail by 



Calvet et al. 



(119831 ). For this purpose we used the CO (3— ^2) 



spectra since it has the best signal-to-noise ratio, but note that using any of the three 
transitions render consistent values within a 10%. In the low velocity range we selected 
—12.5 km s _1 , the ambient cloud velocity, as the limiting value between what we consider 
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blue or red-shifted gas. The derived moment urns are given in Columns 5-7 of Table [3) They 
were computed using the observed radial velocities, and thus correspond to strict lower 
limits. To correct for inclination they should be multiplied by a factor (cos?) -1 where % is 
the angle between the flow and the line of sight. 

The average momentum rates can be estimated from the second moment of the 
CO (3— >2) spectra, from the expression 

P = [H 2 / 12 CO] /i m ^F(J,T ex )J^^ dA , (5) 

lobe rchar 

where F(J,T ex ) is the term that multiplies the integral in the right hand side of Eq. (pQ), 
and we select J = 2 as indicated previously. The opacity t v can be estimated from Eq. 
(J3]) in the optically thin limit. The characteristic radius of the flow, r c h ar , is a sensitive 
parameter that regulates the size of the flow and the time scale (after dividing it by the flow 
velocity) over which the momentum has been deposited in the pre-stellar core. In Eq. (jSJ), 
we consider r char as the projected radius, taken as the length subtended by 30" (0.23 pc at 
1.6 kpc) for the high-velocity flow; and as the length subtended by 20" for the low-velocity 
gas. The momentum rates derived for the outflows, with no correction for inclination, are 
given in cols. 8-10 of Table [3j A correction for inclination will require to multiply by the 
factor (cos i)~ 2 sin i. 

Toward each of the lobes of the SE-NW flow we observe wing emission at both 



blueshifted and redshifted velocities. This a 



using the expression (jCabrit Bertout 



1990 



lows us to make an estimate of its inclination 



19861 ) 



R+l 1 

tarn 



#-ltan0 ' 

where R is the ratio between maximum observed blueshifted and redshifted velocities 
toward a single lobe, and 9 is half of the total aperture angle of the outflow. We estimate 
the later from the 2-MASS _fT s -band image, obtaining a value of ~ 23°. Using this 
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and the observed values of R, we derive an inclination angle for the SE-NW outflow of 
80°. Correction factors for the momentum and momentum rate are then 5.8 and 33, 
respectively. The inclination-corrected momentum of the SE-NW flow is 89 M km s _1 and 
the inclination-corrected momentum rate (calculated with Eq. [5]) is 3.0 x 10~ 2 M km s _1 
yr" 1 . 

One of the principal problems of the previous estimation is that one assumes that the 
momentum has been injected into the molecular gas continuously, despite the molecular 
high velocity gas and the radio images evidence that the emission comes from discrete 
episodes of mass ejection. Therefore, the momentum rates calculated are time-averaged 
over the history of mass ejections from the protostar, that probably includes long periods 
of inactivity. 



Guzman et al. 



fl2010h —corrected 



The momentum flux of the ionized jet reported by 
for 80° inclination — is ~ 4.2 x 10~ 4 M Q km s _1 yr _1 , while that of the associated SE-NW 
molecular flow (3.0 x 1O _2 M km s _1 yr _1 ) is approximately 70 times greater. Also, the 
observed momentum of the flow could not have been deposited in the ambient gas by such 
a jet in less than ~10 5 yr, a time larger than the estimated life of the YSO. A possible 
solution to this inadequacy of the jet as the driving source of the molecular outflow could 
be that the protostellar jet is only partially ionized, being its mass and momentum loads 
greater than the estimated ones. 

4.3. Further tracers of shocks 



In addition to the high-velocity CO gas and its high excitation temperature, further 
evidence for the presence of shocks is provided by the detection of SiO emission. SiO is 
thought to b e generated when strong shocks pass thr ough dense molecular gas disrupting 



dust grains (IHartquist et al. 



198G 



Caselli et al. 



19971 ). Therefore, we expect that SiO 



emission should be associated with the outflowing gas. This suggestion is strongly supported 
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by the almost identical shape of the profile of the high velocity outflowing gas emission 
observed in the HCO + (4—7-3) line and the SiO (8—7-7) profile (see next section and bottom 
panel of Fig. |2J). 

The observed velocity integrated SiO emission is ~ 4.6 K km s _1 . Assuming the SiO 
line is emitted in optically thin conditions with an excitation temperature of 75 K, we 
derive a column density of SiO of 7.8xl0 12 cm -2 . Assuming T ex = 100 K only changes this 
value by 3%. The critical density for the SiO(8— >7) transition is ~ 10 8 cm -3 , considerable 
larger than the ambient density, indicating that thermal collisions with ambient H2 is not 
the main mechanism for excitation to this rotational state. The good correlation between 
the SiO emission and the high velocity component suggests that this mechanism is shock 
excitation of the gas. 



There are three OH masers reported in the literature toward the region (jCaswelllll998 



20041 ) . Figure H] shows their positions with green crosses. The closest to the peak of the SE 



blue-shifted lobe corresponds to a 1720 MHz OH maser, which has an LSR ve 



km s 1 . The 1720 MHz OH masers are usually associated with shocked gas ( IBrogan 



ocity o 



-24.5 



20070 



and the location and velocity of this particular one is con sistent w i th it being pumped by 



the shock produced by the high velocity blueshifted flow. ICaswelll ( 12004] ) reported Zeeman 



splitting on this maser implying a ~ 10 mG magnetic field, the latter probably enhanced 
through shock compression. We also note that the excitation temperature of the outflowing 
gas — derived in §4.2.11— is above the minimum temperat ure required to pump this maser 
collisionally (T w 90 K, 



Elitzur 



1976 



Hoffman et al. 



20031 ). We conclude that this maser is 



associated with the IRAS 16562—3959 core and to the shocked gas in the SE-NW ou tflow, 



and not located farther away at ~ 2.5 kpc as previously suggested by 



Caswell! (J2004J). 



The other two masers correspond to 1665/1667 MHz OH masers: one is coincident with 
the central jet source (Vqh = —12.7km s _1 ) and the other, located about 15.6" NE from 
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the jet (Vqh = —15.5 km s 1 ), is associated with a mid infrared so urce in the field, uns een 



Guzman et al 



20101 ) 



in 2MASS but conspicuous in Spitzer 8.0/im band (see Fig. 4 from 



4.4. Infall of molecular core 



To derive infall parameters we mo deled the HCO + (4— >3) profile using the simple 



analytic model of contracting clouds of 



Myers et al. 



(119961 ) . We first subtracted a broad 



Gaussian profile, shown as a dash-dotted line in the upper panel of Fig. [21 to take into 
account the contribution of the high velocity outflow. We note that the shape of the 
subtracted Gaussian is very similar to the profile of the SiO(8 — > 7) emission which is 
expected to come mainly from the outflowing gas. This is illustrated in the bottom panel 
of Fig. [21 which shows the SiO(8 — > 7) spectrum and superimposed the broad Gaussian 
profile fitted to the HCO + (4— >3) emission (scaled by a factor of 1/14). 

The simple model can very well fit the observed profile (see top panel of Fig. [2]). The 
parameters of the best fitting collapse model are: an infall velocity (T4n) of 0.35 km s -1 , a 
velocity dispersion (a) of 1.54 km s -1 , an optical depth (r ) of 3.8, a kinetic temperature of 
32 K, and a rest-frame velocity of the collapsing envelope of —11.4 km s -1 . This velocity 
is about 1 km s~ l grea ter than the ambient cloud velocity adopted here, but we note that 



Bronfman et al. 



(119961 ) also reported a velocity of —11.6 km s 1 in the CS(2— >1) transition 



detected towards this source. 



In order to estimate the envelope accretion rate we consider the mass dis tribution 



model for the parent cloud derived from dust emission in iGuzman et al. 



(120 10l ). Assuming 



that the extent of the contracting gas along the line of sight is comparable to the size of the 
cloud of about 0.7 pc FWHM (1.5 ' at 1.6 Kpc) and using the derived values of the infall 
speed, molecular density, and core size, we obtain a mass infall rate M in of ~ 3.7 x 1O~ 4 M 
yr " 1 , a value one order of magnitude below the estimation from the SED fitting reported 



by 



Guzman et al 



(]2010l ). The value derived for the infall velocity is also about one order of 



-17- 



magnitude below the free-fall velocity expected for the entire cloud, suggesting a retarded 
collapse. 

5. SUMMARY 

We undertook molecular line observations toward the IRAS 16562—3959 massive and 
dense core in three CO transitions: CO (3— >2), CO (6— >5), and CO (7— >Q), using APEX, with 
angular resolutions ranging from ~8" to ~17". Additionally, HCO + (4— >3), C 18 0(3— >2) and 
SiO(8— t-7) observations were performed towards the center of the core. The main results 
and conclusions presented in this paper are summarized as follows: 

1. High velocity molecular gas was detected in the three CO transitions, spanning a 
range in radial velocity of ~ 60 km s _1 . 

2. The morphology of the high velocity emission, best shown in the CO (6— >5) and 
CO(7->6) lines, is quadrupolar. We conclude that this morphology is due to the 
presence of two collimated bipolar outflows, one lying in the SE-NW direction and 
another in the N-S direction. The more extended bipolar flow is in the SE-NW 
directio n, consistent with be ing excited by the ionized jet detected toward this 



MYSO flGuzman et al 



2010l ). Extended -K^-band emission probably tracing excited 



H 2 -2.12/im is also associated with the SE-NW flow. 

3. From an analysis of the emission in the three CO lines we derive that the excitation 
temperature of the high-velocity gas, assumed optically thin, is high, with values of 
~ 140 in the lobes of the SE-NW outflow and ~ 90 in the lobes of the N-S outflow. 
This suggests that the excitation of the gas is related with the acceleration mechanism, 
namely, shock- induced acceleration. 

4. The total mass in the SE-NW and NS outflows are, respectively, 1.9 and 1.3 M Q . The 
momentum and momentum rates derived for the SE-NW outflow, corrected for an 
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inclination angle of 80°, are 89 M km s _1 yr -1 and 3.0 x 10~ 2 M km s' 1 yr -1 , 
respectively. These values are characteristics of flows driven by young massive stellar 
objects with high luminosities. 

5. The molecular core in which the outflow is embedded presents evidence of being 
in gravitational contraction as shown by the blue asymmetric peak seen in the 
HCO + (4— >3) transition observations. The derived mass infall rate is of the order of 
a few times 1O~ 4 M yr _1 and the infall velocity is ~ 0.35 km s -1 , which is a factor 
of ~ 10 below the velocity expected from a free-fall collapse. This implies that the 
collapse has been retarded and not gravitationally dominated. 

The authors gratefully acknowledge support from CONICYT through projects 
FONDAP No. 15010003 and BASAL PFB-06. This publication is partly based on data 
acquired with the Atacama Pathfinder Experiment (APEX). APEX is a collaboration 
between the Max- Planck- Inst it ut fur Radioastronomie, the European Southern Observatory, 
and the Onsala Space Observatory. This publication makes use of data products from 
the Two Micron All Sky Survey, of the GLIMPSE-Spitzer database, and of the Red MSX 
Source survey database at www.ast.leeds.ac.uk/RMS. 
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Table 1. Observational Parameters 



Line 


Frequency 


Tel. 


Beam 


Vmb 


Map 


Spacing 


Av 


Noise 




[GHz] 




["] 






["] 


[km s _1 ] 


[K] 


CO(3^2) 


345.796 


APEX 


17.3 


0.73 


5x5 


20" 


0.423 


0.1 


CO(6^5) 


691.627 


APEX 


8.8 


0.47 


29 x 21 


4" 


0.318 


0.5 


CO(7^6) 


806.898 


APEX 


7.7 


0.43 


27 x 19 


4" 


0.272 


1.0 


SiO(8^7) 


347.331 


APEX 


17. 


0.73 


3x3 


20" 


0.423 


0.05 


HCO+(4^3) 


356.734 


ASTE 


22 


0.7 


Central Pos. 




0.11 


0.1 


C 18 0(3^2) 


329.331 


ASTE 


23 


0.7 


Central Pos. 




0.43 


0.05 



Table 2. Parameters of High- Velocity and Ambient Gas 



Lobe 


V range 


N CO 


T 

x ex 




(km s- 1 ) 


(10 16 cm" 2 ) 


(K) 


Blue SE 


[-46.0, -22] 


8.72 


145 


Red NW 


[-3, 14.0] 


6.06 


120 


Blue N 


[-46.0, -22] 


5.70 


76 


Red S 


[-3, 14.0] 


3.37 


88 


Ambient 


[-22, -3] 


11.30 


54 
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Table 3. Outflow Parameters a 



Mass Momentum 13 Momentum rate b 

Outflow (M ) (Mekms^ 1 ) (10" 4 M km s" 1 yr" 1 ) 

HV LV Total HV LV Total HV LV Total 



SE-NW outflow 



Blue SE 0.39 0.61 1.00 5.8 2.9 9.7 4.3 1.2 5.5 

Red NW 0.27 0.65 0.92 3.5 3.1 6.6 2.2 1.3 3.5 

Total 0.66 1.26 1.92 9.3 6.0 15.3 6.5 2.5 9.0 

NS outflow 



Blue N 0.25 0.49 0.74 3.4 2.4 5.8 3.3 1.0 4.3 
Red S 0.15 0.43 0.58 1.9 2.0 3.9 1.7 0.8 2.5 

Total 0.40 0.92 1.32 5.3 4.4 9.7 5.0 1.8 6.8 

a HV and LV denotes parameters of the high-velocity and low-velocity flows, 
derived as described in §4.2.21 

b Not corrected for inclination. 
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-40 -20 20 




Offset-RA (") 

Fig. 1.— Spectral grid of the CO (3^2) (red), CO (6^5) (green), and CO (7^6) (blue) 
emission observed towards IRAS 16562—3959. The grid spacing is 20". Offsets are from the 
radio source reference position at a 2 ooo = 16 h 59 m 41.6r!0, 5 2 ooo — — 40°03 / 43". Velocity scale 
ranges from —45 to 25 km s _1 . The dashed line in each spectrum marks the ambient velocity 
adopted of —12.5 km s _1 . 
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Fig. 2. — Spectra observed toward the peak position of IRAS 16562—3959. Top: 
HCO + (4— >3) spectrum. The dashed line shows the best-fit using a collapsing envelope model 
from Myers et al. (1996) plus an outflow component (dot-dashed line). Middle: C 18 0(3— >2) 
spectrum. The dashed-line indicates the adopted ambient cloud velocity of —12.5 km s _1 . 
Bottom: SiO(8— >7) spectrum. The dot-dashed line is not a fit but shows the outflow profile 
fitted to the HCO + (4-)-3) spectrum scaled by 1/14. 
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)(3-2) IRAS 16562-3959 




H ' h 

30(6-5) 
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00(7-6) 




Blue: -46.0 to -22.0 km s 
Red: -3.0 to 14.0 km s" 1 



16 59 45 16 59 40 

Right Ascension (J2000) 

Fig. 3. — Contour maps of the velocity integrated CO line wing emission towards IRAS 
16562—3959. Blue lines represent emission integrated over the velocity range —46 < vi sr < 
—22 km s _1 , which is blueshifted with respect to the ambient velocity of —12.5 km s _1 , and 
red dashed lines emission integrated over the velocity range —3 < vi sr < 14 km s _1 , which 
is redshifted with respect to the ambient velocity. The star marks the position of the jet 
source. Contour levels are 20, 30, 40, 50, 60, 70, 80, and 90% of the peak emission. Top: 
CO(3— >2) emission. Peak blueshifted emission: 49.0 K km s -1 . Peak redshifted emission: 
26.4 K km s" 1 . Middle: CO(6— >5) emission. Peak blueshifted emission: 85.2 K km s _1 . 
Peak redshifted emission: 49.6 K km s _1 . Bottom: CO(7— >6) emission. Peak blueshifted 
emission: 85.1 K km s" 1 . Peak redshifted emission: 54.3 K km s _1 . 
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Fig. 4. — Map of the CO (6— >5) velocity integrated blueshifted (blue contours) and redshifted 
(red contours) wing emission, as in Fig. [31 overlayed with a map of the 8.6 GHz emission 
(black contours) showing the string of radio emission described in Guzman et al. (2010). 
The red and blue dashed-line boxes mark the regions of the sky where the CO emission was 
integrated to determine the parameters of the red and blue lobes, respectively. The green 
crosses marks the position of the OH masers. 
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17 h 00 m 00 s 16 h 59 m 50 s 16 h 59 m 40 s 16 h 59 m 30 s 

Right Ascension (J2000) 

Fig. 5. — Grey scale: i^-band 2MASS emission. Overlaid are contours of the blueshifted 
(blue contours) and redshifted (red contours) CO (6— >5) emission. The approximate direc- 
tions of the SE-NW and NS flows are indicated by dashed lines. 
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Fig. 6. — Velocity integrated CO emission versus upper rotational quantum level J. The 
range of velocity integration for the different features is indicated in Column 2 of Table [2j 
The outflow emission from the lobes has been averaged over the zones indicated in Fig. HI 
Dashed lines represent fits to the data using Eq. (T1J and optically thin conditions. The 
derived column densities and excitation temperatures are shown in Table [21 and the latter 
are also indicated in each panel. Top: Wing emission from the blueshifted lobe of the SE-NW 
outflow. Middle-top: Wing emission from the redshifted lobe of the SE-NW outflow. Middle: 
Wing emission from the blueshifted lobe of the NS outflow. Middle-bottom: Wing emission 
from the redshifted lobe of the NS outflow. Bottom: Ambient cloud emission, measured 
from spectrum observed at offsets ARA = +40", ADec = —40". 



